Co-circulation of infectious diseases on networks 
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We consider multiple diseases spreading in a static Configuration Model network. We make 
standard assumptions that infection transmits from neighbor to neighbor at a disease-specific rate 
and infected individuals recover at a disease-specific rate. Infection by one disease confers immediate 
and permanent immunity to infection by any disease. Under these assumptions, we find a simple, 
low-dimensional ordinary differential equations model which captures the global dynamics of the 
infection. The dynamics depend strongly on initial conditions. Although we motivate this article 
with infectious disease, the model may be adapted to the spread of other infectious agents such 
as competing political beliefs, rumors, or adoption of new technologies if these are influenced by 
contacts. As an example, we demonstrate how to model an infectious disease which can be prevented 
by a behavior change. 



Introduction Networks have captured the attention of 
many scientists. One of the primary interests is in under- 
standing how network structure governs the behavior of 
dynamic processes spreading on networks [6, 16, 32J. This 
is complicated by difficulty in deriving analytic models, 
limiting our understanding of the dynamics. Sometimes 
we can gain insight into long-term behavior without un- 
derstanding dynamics, but in many cases the interme- 
diate dynamics governs the long-term outcome. This is 
particularly significant when competing processes are oc- 
curring in the network. In this article we study the si- 
multaneous spread of two competing diseases in a Con- 
figuration Model network. Although we focus on disease, 
other competing "infectious" processes, such as a change 
in behavior in response to a disease [5] , spread of beliefs 
in a voter model [B], and "viral marketing" of compet- 
ing technologies [2] have been studied, and the approach 
introduced here can be adapted to these applications. 

In this article we derive a low-dimensional system of 
equations capturing the dynamics of competing diseases 
spreading simultaneously in a Configuration Model net- 
work. We apply the model to investigating possible out- 
comes of co-circulating diseases. Prior studies have thor- 
oughly analyzed the effect of network structure such as 
degree distribution [3J [T7J [551 I33J and heterogeneities in 
infectiousness and/or susceptibility [131 HH [201 E3] on 
disease spread. Recent work gives insight into the role 
of partnership duration [2"4Tf2"6l 135] . Other investigations 
focus on the role of clustering HU [El HFJ EH E21 [311 
[331 131], with limited predictive success. 

Models of interacting diseases typically neglect net- 
work structure {e.g., CDS] and many others). Until re- 
cently, models of a single disease spreading through a 
network have relied on approximation [8] or been re- 
stricted to final size calculations under the assumption 
of an asymptotically small initial fraction infected [55] , 
Extending these approaches to competing diseases [101 
IT31 [3D] does not allow us to measure the effect of dy- 
namic interactions, and so results are limited to special 
cases in which these interactions are not important, such 
as when one disease spreads before the other. 

Our method can be easily adapted to more than 2 dis- 



eases and allows for arbitrarily large initial conditions. 
We validate the system by comparison with simulation. 
Using our equations, we are able to identify the sear- 
ings which separate different regimes. We discuss these 
regimes and introduce possible generalizations. 

The basic model We assume that two diseases spread 
in a Configuration Model network [29] (also called a 
Molloy-Reed network [57] ) with degree distribution given 
by P(k). For disease 1 transmission along an edge has 
rate /3i and recovery of infected individuals has rate 71 . 
For disease 2 the rates are /32 and 72 . A node infected by 
either disease gains immunity to any further infection. 

Our approach is similar to that of [T3] and is based 
on [24) . We will focus our attention on a test individ- 
ual u (described more fully in the appendix and (23)1. 
a randomly chosen individual in the population. We as- 
sume that the aggregate population-scale spread of the 
diseases is deterministic. Under these assumptions, the 
probability the test individual has a given infection status 
equals the proportion of the population with that status. 
Thus by calculating the probability a test individual has 
a given status, we immediately know the proportion of 
the population with that status. 

We make one change to the test individual u: we pre- 
vent it from causing infections. This keeps the status of 
its partners independent of one another without affecting 
its own status, and so it has no effect on our calculations 
of the proportion of the population in each state. An 
alternate argument for why this change has no impact 
is that we have assumed the dynamics are deterministic, 
while the timing of when (or even if) u is infected is a ran- 
dom variable. Thus the infections u would cause cannot 
have any macroscopic impact on the disease dynamics, 
and this modification of u has no effect. 

We take t = to to be our "initial time" . In practice this 
may correspond to the time of introduction of a disease if 
enough individuals are initially infected, or it corresponds 
to a later time at which enough infection is present that 
the spread is deterministic. There are some restrictions 
on how the initial infections can be distributed, discussed 
in the appendix. We choose a test individual u randomly 
from the population (it may have any status). We let 
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FIG. 1: (left) flow diagram for the probabilities the test indi- 
vidual has each status, (right) flow diagram for the probabil- 
ities a partner of the test individual has each status. 



v be a random neighbor of a random test individual u 
which had not transmitted to u by to- We define 9(t) to 
be the probability that at time t, v has not transmitted 
to u. The probability u is susceptible at time t is 

S(t)=iPm))=J2 p (k)S(k,t )8(t) k 

k 

where S(k,to) is the probability an individual of degree 
k is susceptible at t = to- We take I\ and I2 [resp R\ and 
R2] to be the probabilities that u is infected with [resp 
has recovered from] the corresponding disease. 

To calculate the change of 9, we must know more about 
the probability v is in any given state. We define 4>s(t) to 
be the probability v is susceptible, 4>i,i(t) and <j>j %{t) to 
be the probabilities that v is infected has not transmitted 
to it, and 4>r to be the probability v is recovered but did 
not transmit (we do not need to distinguish which disease 
infected v). Then 6 — <f>s + 4>i,i + 4>i.2 + 4>r- 

We calculate 4>s similarly to S. If v is initially suscep- 
tible we find the probability it has degree k by counting 
all edges of initially susceptible individuals of degree k: 
NkP(k)S(k,to) and dividing by the number of all edges 
of initially susceptible individuals J2k NkP(k)S(k,to) 
(N is population size). If v has degree k and was initially 
susceptible, the probability v is still susceptible is 9 k ~ 1 
(because u is prevented from transmitting to v). This 
leads to the conclusion that v is susceptible with proba- 

bility M t) = Mto)^§§^ = 

Fig. [l] gives flow diagrams which yield our equations. 
Each box represents a compartment, and arrow labels 
represents probability flux from one compartment to an- 
other. The fluxes from the I compartments to the R 
compartments represent recovery of u. The fluxes from 
the 4>i compartments to the 4>r compartments [resp 1 — 0] 
represent flux due to recovery of v prior to transmitting 
[resp transmission prior to recovery]. The fluxes from S 
and 4>s are found by differentiation of S and 4>$ in time, 
using 9 = —Pi4>i t x—p2 < t>i,2-, and assigning the appropriate 
proportion to the appropriate compartment. 



From the diagram, we find 

9 = -Mi tl - (i) 

ip" (9) 

4>I,m = - (Pm +7m)</>/,m + Pm&I.m&S (*o) 7T77\ (2) 

1p'{l) 

S = (3) 

I m = Pm<f>I,m1p'(8) ~ InJm (4) 

Rm = ^vnJm (5) 

The subscript m takes the values 1 and 2. The single- 
disease, small initial condition limit of these equations 
has been proved exact [5]- These equations capture the 
fact that disconnected components are safe from outside 
introduction. 

Sequential Introduction As an example, we consider 
two diseases spreading in a network of 10 6 individuals 
with Poisson degree distribution of mean 5. For both, 
7 = 1 but Pi = 0.5 and fa = 1. In simulations, we 
introduce disease 1 into 30 random individuals at t = 0, 
and disease 2 into 30 random individuals at t — 5.25. 

Our deterministic equations do not apply while either 
disease has a small number of infections. We take obser- 
vations at to = 4 (after the first, but before the second 
disease) and to — 7 (after both are established) to ini- 
tialize our equations. Fig. [2] compares calculation with 
simulation. Comparing the to — 4 calculation with the 
to = 7 calculation shows the effect of the second disease. 

Simultaneous Introduction We now consider the si- 
multaneous introduction of two diseases and assume that 
the initial numbers infected are large enough that the 
dynamics are deterministic. In our example, we take 
fa = 1.2, 7i =4, fa = 0.2, and 72 = 0.25. Dis- 
ease 1 tends to spread more quickly, but disease 2 has a 
higher probability of transmission prior to recovery. At 
t = 0, we infected a randomly selected proportion of the 
population pi with disease 1, and a proportion p2 with 
disease 2. This gives S(k, 0) = 1 — p\ — p2 for all k, 
h = Pi, h = P2, (f>i,i = Pi, and 0/, 2 = P2, with no re- 
covered individuals. In our population, the degree of each 
node is assigned uniformly from the integers 1 through 
9. We use our equations to calculate the final proportion 
infected by each disease, shown in Fig. [3j 

There are several distinct regimes we can identify in 
Fig. [3] If p 2 is 0(1) and p 1 small, or if p x is 0(1) 
and P2 small, the disease with the large initial condi- 
tion spreads and effectively infects everyone simply be- 
cause a large fraction is initially infected. The other 
disease cannot spread in the "residual network" left 
behind. If neither p\ or p2 is initially large, other 
regimes are seen. To analyze them, we first note that 
when both diseases are small, they grow at exponen- 
tial rates r\ = -(Pi + 71) + fatp" (1) ftp' (1) = 2 and 
T2 = -{P2 + 72) + /W( W(l) = 0.75. 

In the "overlapping epidemic" regime, the slower- 
growing disease 2 begins with a head start. The size of 
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FIG. 2: The spread of two diseases in a population of size 
10 6 with a Poisson degree distribution of mean 5. The first 
disease is introduced with 30 cases at t — 0, and the second 
with 30 cases at t = 5.25. The second strain is more infectious. 
Predictions (dashed) are calculated from observations at t = 4 
(before the second disease's introduction) and t = 7 (shortly 
after). Model and simulation agree well. 



the head start scales so that the two epidemics become 
large at the same time. The value of ln/ 2 — (V2A1) lnli 
is constant during the linear growth phase. For given 
C = In P2 — (r 2 /ri) lnpi, the behavior is universal. 
The diseases grow independently until the linear growth 
phase ends. We can estimate bounds on the regime by 
crudely assuming exponential growth continues forever. 
There is some value of C m i n = In 0.0025 — {ri/r-Cj In 1 = 
In 0.0025 « —6 that corresponds to p% = 1 and p 2 — 
0.0025 which means that the slower growing disease 
would affect less than one percent of the population by 
the time the faster growing disease has fully established 
itself in this approximation. Similarly we take some 
C max = lnl - (r 2 /n)ln0.05 = -(r 2 /n) In 0.05 « 3r 2 /n 
corresponding to p\ = 0.05 and p 2 = 1, which means 
that the slower growing disease will have fully burned 
through the population when the faster growing disease 
is only affecting 5% of the population in this approxi- 
mation. For C < C m in, h becomes large well before 
I 2 . For C > C max , I2 becomes large well before I\. Be- 
tween these values, the epidemics become large at simi- 
lar times and interact dynamically. Further discussion of 
these boundaries is in the appendix. 

There are two "non-overlapping epidemic" regimes. If 
C < C m in, disease 1 becomes large and has an epidemic 
while disease 2 is still exponentially small. If C > C max , 
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FIG. 3: We take a network with the degree of each node 
chosen uniformly from 1 up to 9. A proportion p± of the 
population begins infected with disease I and a proportion 
p2 begins infected with disease 2 with /?i = 1.2, 71 = 4, 
/?2 = 0.2, and 72 = 0.25. Disease I spreads more quickly, but 
has a smaller per-edge transmission probability. The solid 
lines denote the estimated upper and lower bounds of the 
overlapping epidemic regime. 



disease 2 has an epidemic while disease 1 is still expo- 
nentially small. Once one epidemic has finished, it may 
be possible for the remaining disease to spread in the 
residual network. We can derive the threshold condi- 
tion: For simplicity we assume disease 1 spreads first. 
While it spreads, <f>j 2 remains negligible. By looking 
at the relative fluxes out of (f>j t i we conclude that af- 
ter disease 1 has completed its spread but before dis- 
ease 2 is significant, (f>R = 71(1 — 6)/P\. We can as- 
sume <f> s (0) = I and (}> s (t) = ^'{e)li)'{l). Since <f> IA 
and 4>i 2 are effectively zero, we conclude 9 = 4>s + 4>R — 
+ 71 (1 " 0)1 Pi, yielding 



Pi +71 



7i 



no) 

V>'(i) 



Then R± = 1 — S(6). For disease 2 to spread, we re- 
quire 4> Ij2 > 0, which implies P2/W2 + J2) > i)'(l)/tp"(9) 
where 9 comes from the above equation. A similar calcu- 
lation would lead to a final size for i? 2 , so in this case we 
can calculate the final outcomes of the epidemics without 
calculating the dynamics. We note that the threshold we 
have found requires that /3 2 /(/3 2 + 72) be greater than 
Pi /(Pi +71) by a nonzero amount. It is not enough for 
the second disease's transmission probability to simply 
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be larger than the first, it must be well above that of 
the first. This is because even once disease 1 peaks and 
begins to decrease, 9 will continue to decrease further, 
so disease 2 encounters a population that is well below 
the threshold for disease 1. The threshold condition for 
disease 2 to invade has been identified previously: [30] 
derived it under the assumption of a second disease intro- 
duced after the first disease had spread, and [13] derived 
it in the special case that the system was in the non- 
overlapping epidemic regime for which the faster growing 
disease would "win." 

Thus our analysis shows that if two diseases are intro- 
duced at small levels into the network, then the possible 
regimes can be understood by looking at the exponential 
growth rates r\ and r 2 of the diseases. Without loss of 
generality, we can assume r% > ri , the first disease grows 
faster. If disease 2 has a sufficiently large head start, 
there will be a "non-overlapping epidemic" regime: the 
population will experience an epidemic of disease 2 un- 
affected by disease 1. If the infectiousness of disease 1 is 
large enough it will cause its own epidemic after disease 
2 has finished. We can calculate the final size of each epi- 
demic without requiring the full dynamic calculation. If 
we take a smaller head start for disease 2, there is an over- 
lapping epidemic regime in which the two diseases pro- 
duce interacting epidemics. To calculate the dynamics of 
these epidemics, we require the dynamic equations ([!])- 
([5]). The possible final sizes depend on the details of 
the interactions, and there appears to be no simple ex- 
pression for the final size. If the head start for disease 
2 shrinks further, we enter another non-overlapping epi- 
demic regime in which disease 1 has the first epidemic. 
Again, it is possible for disease 2 to later have an epidemic 
if its transmission probability is sufficiently larger than 
that of disease 1. At most one of the non-overlapping 
epidemic regimes can have epidemics for both diseases. 

Generalizations This model may be adapted to other 
"infectious" agents such as the spread of a rumor or a 
new technology through a social network. As an exam- 
ple of its flexibility we consider a disease which can be 
avoided through behavior modification. We assume con- 
tact with an infected individual transmits infection at 
rate /3. However, we allow that if u is in contact with an 
infected individual, then at rate So u changes its behav- 
ior. If it is in contact with an individual who has changed 
behavior then u changes behavior at rate 5b- This leads 
to the flow diagram in Fig. [4] 

In Fig. [8] we show how behavior change modifies epi- 
demic outcomes. If the disease is introduced in a small 
number of individuals and behavior change spreads suffi- 
ciently faster than the disease, then the behavior change 
prevents the disease from having a large-scale epidemic. 

Summary We have introduced an analytic model that 
calculates the simultaneous spread of two infectious dis- 
eases in a configuration model network. Our model 
is low-dimensional regardless of the degree distribution. 
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FIG. 4: The flow diagram for disease spread with behavior 
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FIG. 5: Final sizes of epidemics with behavior change. We 
use networks of the same structure as in Fig. [3] The disease 
spreads with ,9 = 2, 7 = 1. The values of 5 b and 5d are 
varied. For the initial condition, no individuals have changed 
behavior and a proportion 10 -6 is initially infected randomly. 



Using this model, we are able to calculate the effect of 
interactions between diseases in regimes that are inacces- 
sible to analytic theories that do not include dynamics. 

This model can easily be generalized to a range of other 
"infectious" processes. We have shown its application to 
behavior changes in response to a disease. 
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Appendix 

Detailed derivation of equations We begin with a con- 
figuration model network, and consider two infectious 
diseases which infect nodes of the network. We assume 
infection by either disease confers immediate and com- 
plete protection from any future infection by any disease. 
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We index the diseases by 1 and 2. Individuals infected 
by disease i = 1, 2 transmit to their partners at rate /3j, 
and recover at rate ji. 

We assume that at the initial time to enough individu- 
als are infected that the disease spreads deterministically 
(at the aggregated population- level), and that the prob- 
ability an individual of degree k is initially infected is 
S(k,to). We must make an assumption about which in- 
dividuals are initially infected. Namely, if we consider an 
initially susceptible individual u, no information we have 
about u at time to tells us anything about the status of 
its partners. This assumption is satisfied if we initially 
infect a random subset of the population, or even if the 
disease has been spreading for some time before to. This 
assumption is violated if we select high-degree individuals 
and preferentially infect their partners. 

The test individual We now introduce the concept of 
a test individual. This concept is described in more detail 
in [23], and it allows us to simplify our calculations. 

Whether a given individual u is infected at any given 
time is a random variable. However, if we make the as- 
sumption that the disease spreads deterministically at 
the aggregated population-level, then whether or not a 
given individual is infected at any given time cannot have 
any impact on the aggregated scale. If it did, there would 
be stochastic effects visible at the population scale. 

This observation allows us to decouple the status of u 
from the dynamics of the epidemic in the sense that we 
can ignore any feedback from u on the epidemic. To make 
this mathematically rigorous, we simply allow u to be- 
come infected and for its infection to proceed as normal, 
but we disallow any transmission from u to its partners. 
This keeps the status of partners of u independent. 

The probability that u is susceptible equals the prob- 
ability that none of its partners has transmitted to it 
(under the assumption u does not transmit to its part- 
ners) . To calculate the proportion of the population that 
is susceptible, infected, or recovered, we assume that u 
is randomly selected from the population and prevented 
from infecting its partners. We call u a test individual. 
The probability u has a given status equals the propor- 
tion of the population that has that status. 

Deriving the flow diagrams We define 9(t) to be the 
probability that a random neighbor of u which had not 
transmitted to u by time t = to still has not transmitted 
by time t. Then if u has degree fc, the probability it 
was initially susceptible is S(k,t ), and the probability 
it is still susceptible is S(k, to)0(t) k . Averaging over all 
possible values of A:, we have 

S(t) = 1>(9(t)) = P(k)S(k,toW) k 

k 

The value of S reduces over time as infections occur. 
Mathematically this appears as a reduction in 9. We 
must calculate how quickly 9 changes, and how much of 
that change is due to each disease. This will allow us to 
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FIG. 6: Flow diagram that leads to the evolution of 9. A label 
along an edge gives the flux of probability along that edge. 



calculate how much of the reduction in S should go into 
each disease's infected class. 

We now look for the change in 9. We define v to be a 
random neighbor of u which had not transmitted to u by 
to- Then 9 is the probability v has not transmitted to u 
by time t. As our initial condition, we have 6 (to) — 1. We 
divide 9 into four compartments. We take <f>s to be the 
probability that v is still susceptible, tfii.i the probability 
v is infected with disease 1 but has not transmitted to 
u, 4>i,2 the probability v is infected with disease 2 but 
has not transmitted to u, and <pn the probability that v 
is recovered (from either disease) and did not transmit 
during infection. These add up to 9: 9 — 4>s + 4>i,i + 
4>i.2 + 4>R- The probability that v has transmitted to u 
is 1 — 9. The flow between these compartments is shown 
in figure [6] 

We first find the rate of change of 9. It is relatively 
straightforward to see that 
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This is because the only path for 9 to decrease is through 
v transmitting to u, which requires that v be infected 
(and not yet transmitted to u). 

At to, we have 4>s(to) is the probability a random 
neighbor of u is still susceptible (given that it has not 
transmitted to u). We take this as an input value. The 
value of 4>i,i, 4>i,2i and 4>r are similarly all input from the 
conditions at to- We can explicitly calculate the value of 
cf>s at later times, if we know 9. To do this, we find the 
probability distribution for degree of v, and then calcu- 
late the probability that no neighbor of v has transmitted 
to v. 

At time to, the edge joining u to v is simply an edge 
from u to a random neighbor that is susceptible. The 
probability that edge connects to a degree k individual is 
proportional to the number of edges that all susceptible 
individuals of degree k have, NkP(k)S(k,to). The nor- 
malization factor is the total number of all edges of sus- 
ceptible individuals, J2k< Nk' P(k , )S(k / ,t ). The proba- 
bility that v is still susceptible at time t is 9(t) k ~ 1 . So the 
probability of having a degree k susceptible neighbor at 
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time t is NkP(k)S(k,t )6(t) k - 1 /£ fc , Nk' P(k')S{k' ,t ). 
Summing over all possible k, and cancelling N, we ar- 
rive at Ek kp ( k ) s ( k ^oW(t) k ^/Ek' k'P(k')S(k',t ) = 
ip'(9(t))/i(j'(l). So given that v is initially susceptible, 
the probability that v is susceptible at a later time t is 
ip'(9(t))/ip'(l). Since the probability v is initially suscep- 
tible is 4>s (to), we conclude 
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FIG. 7: Flow diagram leading to equations for the proportion 
of the population in each compartment. 



The rate of change of (ps is simply 0s(i o )0V , "($)/'0'(l)- It 
is straightforward to see that the amount that goes from 
0s to 0/,! is 0s(io)/?i07,iV'"(0)/'0'(l) an d the amount 
going into 07 )2 is /3 2 0/, 2 <0"(0)/V ''(1). 

So in figure [6j we have expressions for the flux along 
each edge except the edges into 0r. These edges are 
straightforward, because the recovery rate for disease 1 
is 71 and the recovery rate for disease 2 is 72. So the 
total flux from (f>i t i to <pR is 7i07,i and the flux from 0/ 2 

is 7207,2- 

Using the flows in figure [6j we can arrive at a coupled 
system for 9, 0/ 1 and 0/ 2- It is 
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^(1) 

(Jh + 72)0/,2 + /320/,20s(io)^77TY 

(1) 



with 0(i o ) = 1 and 0s(i o ), 0/,i(io) 
the initial state of the population. 



and 0/,2 (to) given by 



These equations govern the spread of the disease 
through the network. However, they are not the usual 
variables of interest. Typically we want to know the 
proportion susceptible, infected, or recovered. Figure [7] 
shows a flow diagram governing the proportion of the 
population in each state. We can use this to recover S, 
Ji, I2, R%, and i? 2 - As we noted above, S(t) = ip{6). 
The flux into I\ can be calculated to be /3 1 07. 1 -0'(0), and 
the flux into I2 is /3207,2V , '($)- The fluxes from each of 
these into the recovered states are 71/1 and 72/2- We 
distinguish the two recovered states because we will fre- 
quently be interested in the total proportion infected by 
each disease. We could have similarly subdivided 4>r into 
two compartments, but it would not provide any infor- 
mation that is useful here. 



Thus our final system of equations 
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Regime Analysis There are several regimes that can 
be identified. When the cumulative number of infections 
is small enough that S and 0s are approximately 1, then 
we can neglect nonlinear terms. Which regime is ob- 
served is determined by the relative sizes of the two epi- 
demics when the linear approximation breaks down. We 
focus our attention on regimes for which the linear ap- 
proximation is valid at the initial time. 

When the initial condition is small, the linear terms 
dominate and the epidemics grow (or decay) exponen- 
tially. The equations governing the epidemics are un- 
coupled in this regime. Physically, this means that com- 
petition for nodes is so weak that it can be neglected. 
In fact a stronger statement is true: in addition to not 
having inter-disease competition, there is no intra-disease 
competition: a transmission path is very unlikely to en- 
counter a node infected along another transmission path. 

Assume that the exponential growth rates of the two 
diseases are r\ and r 2 , with n > r 2 . This continues un- 
til one of them becomes large enough that competition 
begins to appear. At this point the growth of both epi- 
demics starts to slow. If the difference in epidemic sizes 
is large enough at this point, then this first epidemic 
will not be slowed by the much smaller other epidemic. 
The dynamics will proceed as if there were just one dis- 
ease spreading. Eventually the susceptible population 
will decrease, the disease will peak and eventually decay 
away, all with the second disease negligible. Once this 
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decay has occurred, there will be a "residual" network. 
The second disease will continue to spread along this net- 
work. If the residual network is well-enough connected 
(and the second disease sufficiently infectious), the sec- 
ond disease can continue to grow, and then it experiences 
its own epidemic. 

If the diseases are close enough in size when nonlinear 
terms become important, then both diseases contribute 
a non-negligible amount to reduction in S and 4>s at the 
same time. Thus they interact dynamically. Each dis- 
ease contributes in a non-negligible way to hindering the 
spread of the other. To determine whether this can hap- 
pen, we use a simple balance based on the initial sizes 
and the early growth rates. Typically we might expect 
that one disease becomes large while the other is still 
exponentially small. 

If one disease is sufficiently small when the other dis- 
ease becomes large, then the larger disease will spread 
and cause an epidemic that is effectively the same size as 
it would be in the absence of the smaller disease. Using 
the initial sizes and growth rates, it is straightforward to 
calculate the sizes of the two diseases once nonlinearities 
begin to be significant. For the two diseases to not inter- 
act, the "smaller" disease must remain negligibly small 
until the "larger" disease has finished its epidemic. We 
derive slightly different thresholds for the case where the 
smaller disease is the fast-growing disease or the slow- 
growing disease. To derive the threshold condition, we 
make a crude assumption that the two diseases continue 
spreading according to the linear growth rate. In the 
non-overlapping regimes, the smaller disease must remain 
small throughout the spread of the larger disease. We 
crudely choose to apply our conditions when the expo- 
nential growth implies that the larger disease would have 
reached size 1. We look at the size of the smaller disease 
(assuming exponential growth) at this time. For a given 
observed size for the smaller disease, we reach different 
conclusions if it is the fast or the slow growing disease. If 
it is the fast-growing disease, and we observe 0.05, then 
that means that at previous times it was much smaller, 
and so we would not expect to observe any impact. On 
the other hand if it is the slower-growing disease, and we 
observe 0.05, that means for much of the spread of the 
larger disease it was at about that size, and so we would 
expect to observe some impact. So if the fast- growing dis- 
ease is the small disease, we allow it to be as large as 0.05 
when the slow disease would reach 1. If the slow-growing 
disease is the small disease, we require that it be much 
smaller, choosing 0.0025 = 0.05 2 instead. Our choice of 
threshold is somewhat arbitrary, and influenced by the 
fact that In 0.05 « —3 giving a simple expression for our 
thresholds. Taking this, we can derive the C m i n ~ —6 
and C max « 3r2/ri thresholds in the text. 

Taking pi and p 2 to be the proportions initially in- 
fected with each disease (at random). So long as C = 
In Pi — (fa/rt) In Pi is n ot between —6 and 3r2/rl, then 
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FIG. 8: The flow diagrams underlying the spread of an infec- 
tious disease in the presence of a behavior change. We assume 
behavior change can be triggered by contact with an infected 
individual or with an individual who has already adopted the 
change. We assume behavior change provides complete im- 
munity to disease. 



the two diseases will not interact dynamically. One dis- 
ease will become large and run through its entire epi- 
demic prior to the other disease (possibly) having its 
own epidemic. This defines the "non-overlapping" epi- 
demic regimes. If C lies within this range, then we have 
the overlapping epidemic regime. 



Impact of Behavior change We now consider the 
spread of a single disese through a network, which can 
be prevented with a behavior change. We assume that 
the behavior change gives complete protection from in- 
fection. Once an individual has adopted the behavior 
change, the change is permanent. An individual who has 
changed behavior will transmit that behavior change to 
partners. 



The disease transmits at rate j3 and recovery occurs 
at rate 7. Contact with an infected individual causes 
behavior change at rate 6d- Contact with an individual 
whose behavior has changed causes behavior change at 
rate 8 b ■ The flow diagrams are shown in figure [8] The 
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resulting equations are 

= -{f3 D + 5 D )4> D - 5 b 4>b 

<i>B = (Sb<Pb + 5 D (j) D )(j)s{to) ^ l ^} 

<j>D = l3D(t>D4>s{h) ^j-j - {Pd + S D + ^))4>D 

1 = PdM'{0) - 1 i 

B = {5 B <t>B + 5D<t>DW{0) 

R = «fl 

S = 

The early growth of cf>]j is exponential with rate 
pD<j>s(toW{l)/i>'{l)-(pD + S D +l). ItS D is sufficiently 
large, the infection decays. This corresponds to a balance 
between transmitting disease prior to either recovering 
or transmitting behavior change. When the transmission 
probability is small enough a typical infected individual 
causes fewer than 1 new infection and the disease must 
die out. 

Changing 6b does not alter the early growth rate of the 
disease. So we might anticipate that the disease will be 
able to cause a large scale epidemic. However, if 5b > 0, 
the behavior change itself also leads to an "epidemic". 
If the behavior growth rate is sufficiently large, its "epi- 
demic" occurs while the disease epidemic is still expo- 
nentially small. In this limit, the behavior change will 
dominate the populaion, and the disease remains expo- 
nentially small. We do not see the complementary case 
in which the disease becomes large while the behavior 
change is exponentially small, because as disease inci- 
dence increases it directly induces behavior changes. So 
either the behavior change has its "epidemic" first, or the 
two have overlapping "epidemics" . 
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